NOLTR  69-52 


O 

00 

CO 

<7> 

•sO 


A  TIME-DEPENDENT  ANALYSIS  FOR  QUASI- 
ONE-DIMENSIONAL  NOZZLE  FLOWS  WITH 
VIBRATIONAL  AND  CHEMICAL 
NONEQUILIBRIUM 


By 

John  D.  Anderson,  Jr. 


NOL" 


OCT  6 1963  l 


H 


UNITED  STATES  NAVAL  ORDNANCE  LABORATORY,  WHITE  OAK,  MARYLAND 


CN 

•O 


o 

oc 


O 

z 


ATTENTION 


This  document  has  been  approved  (or 
public  release  and  sale,  its  distribution 
is  unlimited. 


NOLTR  69-52 


A  TIME-DEPENDENT  ANALYSIS  FOR 
QUASI-ONE-DIMENS IONAL  NOZZLE  FLOWS  WITH 
VIBRATIONAL  AND  CHEMICAL  NONEQUILIBRIUM 


Presented  by: 
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ABSTRACT:  A  new  technique  is  presented  for  the  numerical  solution 
of  quasi-one-dimensional ,  vibrational  and  chemical  nonequilibrium 
nozzle  flows  including  nonequilibrium  conditions  both  upstream 
and  downstream  of  the  throat.  This  new  technique  is  a  time- 
dependent  analysis  which  entails  the  explicite  finite-difference 
solution  of  the  quasi-one-dimensional  unsteady  flow  equations  in 
steps  of  time,  starting  with  assumed  initial  distributions 
throughout  the  nozzle.  The  steady-state  solution  is  approached 
at  large  values  of  time.  A  virtue  of  the  present  time-dependent 
analysis  is  its  simplicity,  which  prevails  from  its  initial  physical 
formulation  to  the  successful  receipt  of  numerical  results.  Also, 
the  present  solution  yields  the  transient  as  well  as  the  steady- 
state  nonequilibrium  nozzle  flows.  To  exemplify  the  present  analysis, 
results  are  given  for  several  cases  of  vibrational  and  chemical 
nonequilibrium  expansions  through  nozzles. 
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A  Time-Dependent  Analysis  for  Quasi-One-Dimensional  Nozzle  Flows 
With  Vibrational  and  Chemical  Nonequilibrium 


This  report  presents  a  new  technique  for  the  numerical  solution  of 
quasi-one-dimensional  nonequilibrium  nozzle  flows.  The  governing 
equations  and  numerical  approach  are  discussed  in  detail,  and 
results  are  presented  for  the  cases  of  a  calorically  perfect  gas, 
vibrational  nonequilibrium  and  chemical  nonequilibrium. 
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NOMENCLATURE 

A  =  cross-sectional  area,  also  symbol  for  atomic  species 
A '  =  A/A* 

A*  =  area  of  throat 

aQ  =  reservoir  frozen  speed  of  sound 

a  =  frozen  speed  of  sound,  also  a  constant  in  expression  for  t 

a'  =  a/a0 

cv  =  specific  heat  per  unit  mass  of  mixture 

Cv  =  specific  heat  per  mole  of  species  i 

D  -  dissociation  energy  per  molecule 

E^  =  internal  energy  per  mole  of  species  i,  including  the 
heat  of  formation 

V  *  VEo 

Eq  *  internal  energy  of  the  reservoir  mixture 
ev^b  =  vibrational  internal  energy  per  unit  mass 

evib  *  evib^RTo 

h  =  Planck's  constant 

k  *  Boltzmann  constant 

k_^  =  forward  reaction  rate  constant 

F 

kR^  =  reverse  reaction  rate  constant;  *  Kg 

Ke  =  equilibrium  constant 

L  =  characteristic  length  taken  equal  to  length  of  nozzle 
TTl  =  molecular  weight  of  mixture 
p  =  pressure 

q  =  general  nonequilibrium  variable  in  equation  (5) 

R  =  specific  gas  constant  of  the  mixture 
A  *  universal  gas  constant 
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Subscripts 

i  =  chemical  species  i 

o  =  reservoir  conditions 

A  =  atomic  species 

-  diatomic  species 
Superscript 

j  *  collison  partner  in  equation  (14) 
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INTRODUCTION 

Because  of  the  practical  importance  of  high-temperature  flows 
through  rocket  nozzles  and  high-enthalpy  aerodynamic  testing 
facilities,  intensive  efforts  have  been  made  during  the  past 
decade  to  obtain  relatively  exact  numerical  solutions  for  the 
quasi-one-dimensional  expansion  of  a  high-temperature  gas  through 
a  convergent-divergent  nozzle  when  vibrational  and/or  chemical 
nonequilibrium  conditions  prevail  within  the  gas  both  upstream 
and  downstream  of  the  throat.  (See  for  example,  Refs.  1-6,  10-13, 
23-25.)  All  of  these  nonequilibrium  solutions  involve  steady- 
state  analyses,  and  are  by  no  means  trivial;  however,  the 
evolution  of  such  efforts  has  produced  adequate  and  sophisticated 
techniques  for  the  analysis  of  nonequilibrium  nozzle  flows. 

An  authoratative  discussion  of  these  steady-flow  techniques  can 
be  found  in  Ref.  1. 

The  purpose  of  the  present  paper  is  to  present  a  new, 
alternative  approach  for  the  numerical  solution  of  quasi-one- 
dimensional  nonequilibrium  nozzle  flows.  This  new  technique  is 
a  time-dependent  analysis  which  entails  the  finite-difference 
solution  of  the  quasi-one-dimensional  unsteady  equations  of  change 
in  steps  of  time.  For  specified  equilibrium  reservoir  conditions 
and  a  fixed-nozzle  shape,  the  physical  gas-dynamic  properties 
are  obtained  in  steps  of  time,  starting  with  assumed  distributions 
throughout  the  oozzle.  The  steady-state  solution  is  approached 
at  large  values  of  time. 
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The  present  time-dependent  analysis  considers  vibrational 
and  chemical  nonequilibrium  conditions  both  upstream  and  downstream 
of  the  nozzle  throat;  the  proper  steady-state  critical  mass  flow 
is  automatically  approached  at  large  values  of  time.  In  addition, 
no  oscillations  and  instabilities  occur  in  regions  of  near 
equilibrium  flow,  no  special  procedures  are  required  to  start  the 
solutions  from  equilibrium  reservoir  conditions,  and  very  large 
spacings  between  grid  points  can  be  employed  throughout  the 
entire  nozzle.  (Accurate  solutions  have  been  obtained  for 
convergent-divergent  nozzles  with  area  ratios  of  10  using  as  few  as 
15  grid  points  beginning  at  the  reservoir.)  Also,  the  present 
method  readily  handles  simultaneous  rate  processes  involving  very 
slow  and  very  fast  reactions. 

The  main  virtue  of  the  present  time-dependent  technique  is 
its  simplicity.  The  governing  conservation  equations  are 
directly  solved  by  a  simple,  explicit  finite-difference  scheme. 

The  technique  requires  no  additional  mathematical  methods  to 
overcome  special  difficulties  that  can  occur  in  the  analysis 
of  nonequilibrium  flows.  (See  Ref.  1  for  a  detailed  discussion 
of  these  difficulties.)  Also,  the  present  technique  lends  itself 
to  particularly  simple  programming  for  a  digital  computer. 
Consequently,  the  simplicity  of  the  present  technique  prevails 
from  its  initial  physical  formulation  to  the  successful  receipt  of 
numerical  results. 

A  second  virtue  is  that  the  present  technique  yields  the 
transient  as  well  as  the  steady-state  nonequilibrium  nozzle  flows. 
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A  transient  solution  is  of  interest  in  its  own  right;  for  example, 
the  physical  times  required  for  the  fluid  dynamic  and  chemical 
variables  to  approach  their  steady-state  values  can  be  assessed. 

ANALYSIS 

The  physical  problem  treated  in  the  present  paper  is  that 
of  the  quasi-one-dimensional  motion  of  a  high-temperature  gas 
expanding  from  equilibrium-reservoir  conditions  through  a 
convergent-divergent  nozzle,  where  the  expansion  is  rapid  enough 
such  that  vibrational  and  chemical  nonequilibrium  prevail  locally 
in  both  the  subsonic  and  supersonic  portions  of  the  flow.  (The 
term  "quasi-one-dimensional"  will  be  discussed  in  a  subsequent 
paragraph.)  Vibrational  and  chemical  nonequilibrium  are  the  only 
dissipation  mechanisms  assumed  in  the  flow;  the  effects  of  thermal 
conduction,  diffusion  and  viscous  dissipation  are  assumed  to  be 
negligibly  small.  The  equilibrium  reservoir  conditions  and 
nozzle  shape  are  specified,  and  a  solution  is  sought  for  the 
steady-flow  distribution  of  p,  u,  T,  p,  ev^b  and  chemical  composi¬ 
tion  in  the  x  direction  along  the  nozzle  (see  Fig.  1) .  The  unique 
aspect  of  the  present  analysis  is  that  a  time-dependent  approach 
is  used  to  obtain  the  steady-state  nonequilibrium  nozzle  flow 
variables . 

The  crux  of  the  present  analysis  is  as  follows:  For  specified 
equilibrium  reservoir  conditions  and  a  fixed-nozzle  shape,  initial 
values  of  the  flow-field  variables  p,  u,  T,  e  ^  and  chemical 
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composition  are  assumed  at  equally  spaced  grid  points  along  the 
x-axis,  as  shown  in  Figure  1.  If  all  flow-field  variables  are 
known  at  time  t,  then  at  each  interior  grid  point  new  values  can 
be  obtained  at  time  (t  +  At)  from  the  first  three  terms  of  a  Taylor's 
series  expansion  in  time, 


g(t  +  At)  =  g(t)  +  (||)_  At  +  (^-4) 


(At) 


(1) 


3t‘ 


where  g  signifies  In  p,  In  u,  In  T,  ln(e^^)  ,  and  In  n^,  and  At 

is  a  small  increment  in  time  chosen  to  satisfy  certain  stability 

criteria  discussed  in  a  subsequent  section.  (For  the  present 

numerical  computations,  the  natural  logarithms  of  the  nondimensional 

flow-field  variablec  are  employed  as  the  dependent  variables  in 

order  to  improve  numerical  stability,  and  a  nondimensional  time, 

t'  ,  is  employed  in  lieu  of  t  in  equation  (1)  .)  Starting  with  the 

initially  assumed  gas-dynamic  variables  at  t  =  0,  the  flow  field  is 

subsequently  obtained  in  steps  of  time  from  equation  (1) .  At 

large  values  of  time  (after  many  time  steps,  usually  on  the  order 

of  700  or  more)  the  steady-state  flow  field  is  obtained,  where 
2 

■>4  and  both  approach  zero.  For  the  present  investigation 

this  steady-state  solution  is  the  desired  result;  indeed,  the 

purpose  of  the  present  report  is  to  show  that  the  time-dependent 

approach  is  simply  an  advantageous  means  to  the  end. 

The  time  derivatives  (4.)  and  (•— §■)  which  appear  in 

dt  t  t 

equation  (1)  are  obtained  from  ;\e  unsteady  quasi-one-dimensional 
conservation  equations  (the  independent  variables  are  x  and  t) . 
These  equations  can  be  derived  from  the  consideration  of  a  fixed 
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control  volume  in  the  flow;  they  do  net  necessarily  follow  from  the 
general  three-dimensional  partial-differential  equations  of  change 
(formulated  in  any  good  text  in  gas  dynamics,  such  as  Ref.  14) 
because  of  the  physical  inconsistency  between  the  two  assumptions 
of  one-dimensional  flow  and  variable  cross-sectional  area  (hence 
the  term  "quasi-one-dimensional" ) .  Nevertheless,  the  quasi-one- 
dimensional  equations  are  commonly  used  as  a  satisfactory  approxima¬ 
tion  for  the  variation  of  flow  properties  through  a  nozzle;  this 
approximation  becomes  more  accurate  as  the  variation  of  A  with  x 


becomes  sma 

Her. 

For 

a  fixed 

A  =  A 

(x)  , 

the  unsteady  quasi-one- 

dimensional 

equations 

are : 

Continuity : 

a2p 

d  t 

+ 

3 (puA) 
3x 

=  0 

(2) 

Momentum: 

9u 

P3t 

+ 

3u 

PU3*  = 

-  |E 

3x 

(3) 

Energy : 

n3e 

P3t 

+ 

3e 

PU33T  = 

3u 

-  pu 

3  (In  A) 

3x 

(4) 

Rate : 

la 

at 

+ 

u|a  = 

3x 

w(p  ,  T 

f  q) 

(5) 

State : 

p 

52 

pRT 

(6) 

where  q  denotes  a  nonequilibrium  quantity  such  as  e  ^  or  chemical 

composition,  and  w  is  a  function  which  depends  on  the  type  of  rate 

process  under  consideration.  Defining  the  following  nondimensional 

variables,  p'  =  p/pQ,  u'  =  u/aQ,  T‘  =  T/Tc,  x'  =  x/L,  t'  =  t/(L/aQ), 

A'  =  A/A*,  Z  =  In  p',  V  =  In  u',  $  =  In  T' ,  where  the  subscript 

zero  denotes  reservoir  conditions,  L  is  the  length  of  the  nozzle, 

1/2 

A*  is  the  nozzle  throat  area,  and  aQ  =  (y0RqT0)  is  the  reservoir 
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frozen  speed  of  sound,  and  after  replacing  p  in  equation  (3)  with 
equation  (6)  ,  equations  (2)  and  (3)  become 


3V  .  T' 

JT  ~  u 


3  (In  A') 
3x ' 

+  IX,  + 

3X1 

12.1 

TP] 

(7) 

\  /  3<t> 

+  12.  j 

TP ' 

-  u* 

3  V 

dX  1 

(8) 

(Note  that  equation  (8)  applies  to  a  nonreacting  gas,  where  R  = 

constant;  for  a  reacting  gas,  where  R  is  variable,  equation  (8) 

is  slightly  modified,  as  shown  in  a  subsequent  section.) 

In  equations  (7)  and  (8) ,  all  x-derivatives  are  calculated  by 

central  finite  differences  obtained  from  the  known  flow  field  at 

time  t.  Consequently,  equations  (7)  and  (8)  directly  yield  (|^.) 

32a  t 

required  in  equation  (1)  ,  where  g  =  Z  and  V;  ( — *■)  can  be  obtained 

3t2t 

by  differentiating  the  same  equations  with  respect  to  time. 


,  ,  3ZV  .  3  o  ,  .  ,3V  *  , 3 Z  , 

u  i3x, 3t,  +  3X ' 3 t 1 J  'Tt1' 


f  T»  .  ,3%  .  92Z 

'TIP"’  1 3x '  3 1 1  SxT 


sPT tT  +  (H'  +  (!$ 


3y\  i 

W)] 


, ,  3  V  .  3V  3VX 
u  3x 1  3 1 1  +  TP  W} 


However,  the  above  additional  time  derivatives  introduce  cross 

32_ 

derivatives,  which  can  be  obtained  by  differentiation  of 

equations  (7)  and  (8)  with  respect  to  x'. 


3xr5tT 


-«.,2iliap  + Jfv  +^*,  +  (|v( 
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32V 
3  Vat' 


/  T '  i  r  32(t  ,  a2z 

'  v  Trr)  l - o  +  — 


V1 


3x 


,2 


3x 


,2 


/  5 6  .  3Z  v  /3<J> 

^av  axr;  lajT 


av 

ax1 


•)] 


uMj!v  av  2 
ax'2  ^ 


(8b) 


To  complete  the  formulation  of  the  problem,  the  nondimensional 
counterparts  of  equations  (4)  and  (5)  must  be  obtained.  These 
counterparts  have  a  distinct  form  for  each  of  the  three  cases 
treated  in  the  present  paper,  namely  the  cases  of  a  calorically 
perfect  gas,  vibrational  nonequilibrium,  and  chemical  nonequilibrium. 
These  forms  are  presented  in  the  following  sections. 

Calorically  Perfect  Gas 

By  treating  the  case  of  a  calorically  perfect,  nonreacting 
gas  (constant  y) ,  results  have  been  obtained  which  illustrate 
the  purely  fluid  dynamic  behavior  of  the  present  time-dependent 
analysis.  For  a  complete  formulation  of  this  case,  only  the  energy 
equation  needs  to  be  considered  in  addition  to  equations  (7)  and 
(8).  With  e  =  cvT  ar.d  cy  =  R/(y  ~  1)  ,  the  nondimensional  form  of 


equation  (4)  is 


34> 

7tT 


-  (y  -  l)  u' 


av 

ax1 


-  u'  |£r  -  (Y  -  1)  U1 


3  In  A' 

ax' 


(9) 


Equation  (9)  directly  yields  (l^)  required  in  equation  (1) ,  where 

a2q  t  a2_ 

g  =  4>;  as  before,  and  the  resulting  cross  derivative 

are  obtained  by  differentiating  equation  (9)  with  respect  to  t' 

and  x'  respectively. 
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3^4> 


2 

3  V 


,3V  x  ,3V 


3 1 


3  u'>3T^F-+  <Sr)(„T>  + 


3V  3  (In  A’) 
3t 1  cTx"1 


] 


-  r  92  *  +  /3*  \  (W  n 

3x '  3 1 '  (9xr)  (TTr)J 


(9a) 


3  2  <{> 

dX  '  3 1 1 


=  -(Y  -  1) 


3x 


,2 


,3V  x2  .  3 2  (In  A')  .  ,3Vv  3  (In  A')  , 

(Wr)  +  “  +  {W]  3x 1 


-  u'  [ 


324> 


3x 


fl  +  (Itt)  (^t)  1 


f34> 


‘  3x 


(9b) 


Vibrational  Nonequilibrium 

The  present  time-dependent  analysis  has  been  applied  to  the 

vibrational  nonequilibrium  expansion  of  a  pure  diatomic  gas. 

These  conditions  are  illustrative  of  many  practical  nozzle 

applications  where  the  reservoir  temperature  is  high  enough  to 

excite  the  vibrational  energy  mode,  but  not  high  enough  to  result 

in  noticeable  dissociation  of  the  gas.  In  addition,  the  vibrational 

mode  is  assumed  to  be  in  equilibrium  within  itself  but  not  in 

equilibrium  with  the  translational  and  rotational  modes;  i.e., 

a  Boltzmann  distribution  is  assumed  to  exist  within  the  vibrational 

mode.  Consequently,  a  vibrational  temperature  Tv^^,  not  necessarily 

equal  to  the  translational  temperature  T,  can  be  defined.  With 
5 

e  =  Tj-  RT  +  equation  (4)  becomes 


3T 

at 


2 

5 
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3u 

3x 
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3  (In  A) 
3x 


.  I  (— vjb  +  u  !fvib 
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3T 

3x 


(10) 


Defining  =  evit/RT0  and  6  =  ln^evib^'  the  nondimensional  form 

of  equation  (10)  is 
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In  the  rate  equation,  equation  (5),  q  =  and  w  is  assumed  to 
be  the  relaxation  equation  for  a  harmonic  oscillator.  Therefore, 
equation  (5)  becomes 


9e  . . 
vib 

9t 


x  [(evib)eq" 


9e 


evib] 
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vib 


9x 


(12) 


where  t  is  the  vibrational  relaxation  time,  x  =  f(p,T),  and 

(ev^.  )  is  the  equilibrium  vibrational  energy  evaluated  at  the 
eq 

local  translational  temperature  of  the  gas,  (e  . .  )  =  (hv/k)R/ 

viD  eq 

[exp  (hv/kT) -1]  .  Defining  x'  =  x/(L/a0),  the  nondimensional  form 
of  equation  (12)  is 
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Equations  (11)  and  (13)  directly  yield  (42.)  for  use  in  equation  (1)  , 

dt  t 

2  2 

where  g  =  0  and  8.  As  usual,  values  for  (“£2")  and  a^ft  are 
obtained  by  differentiation  of  equations  (11)  and  (13)  with  respect 
to  t'  and  x'  respectively 
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3|evib>ea  _  (hv,  (el-vAT)!^  T, 
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If  t  is  assumed  to  vary  as  ip  =  ae  ,  then 
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Finite-difference  expressions  for  ttt  and  —  can  also  be 
used,  as  discussed  in  a  subsequent  section. 

Initial  values  at  t  =  0  are  assumed  for  p,  u  and  T.  For 
convenience,  the  initial  values  for  e^^  in  the  subsonic  and  throat 
regions  are  equilibrium  values  a.t  the  local  assumed  T,  whereas  at 
some  arbitrary  point  downstream  of  the  throat  ev^  is  usually  set 
equal  to  a  constant  for  the  remainder  of  the  nozzle.  Consequently, 
the  initial  evib  distribution  is  at  least  qualitatively  similar 
to  the  actual  steady-state  nonequilibrium  distribution. 

Chemical  Nonequilibrium 

The  present  time-dependent  analysis  has  also  been  applied  to 
the  case  of  a  nonequilibrium  dissociating  symmetrical  diatomic  gas 

k  j 
kf 

A2  +  M.  t  2A  +  M.  (14) 

3  k  j  3 

R 

where  ML  is  a  collision  partner  (catalytic  body)  ;  in  the  present 
application,  =  A2  and  M2  =  A.  Also,  chemical  nonequilibrium 
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is  the  only  rate  process  considered;  the  vibrational  energy  is 
assumed  to  be  in  equilibrium  with  translation  and  rotation  at  the 
local  gas  temperature.  Hence,  for  simplicity,  no  attempt  is 
made  to  account  for  coupling  between  vibration  and  dissociation, 
even  though  such  coupling  effects  may  be  important  in  some 
practical  applications,^  Letting  subscript  i  denote  a  given 

15 

chemical  species,  e  =  E^,  where  is  the  mole-mass  ratio 

(moles  of  species  i  per  unit  mass  of  mixture)  ,  is  the  molar 
internal  energy  of  species  i  including  its  heat  of  formation,  and 
dE^  =  Cv  dT.  Also,  the  local  frozen  specific  heat  per  unit  mass 

of  the  mixture  is  c  =  E .  n .  C  ,  and  the  local  value  of  the 

vf  1  'i  v^ 

specific  gas  constant  is  R  =  RI,  ri..  In  terms  of  the  above  quantities, 
equations  (3)  and  (4)  become 


8U-  -aU--CpRi|  +  Tp||+T.  ft) 


p  3t  +  pU  3^ 


(15) 
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(16) 


and  equation  (5)  becomes 


-  u 


9x 


(17) 


where  w^  is  the  rate  of  production  of  species  i  by  internal  chemical 

reactions.  Defining  =  r'i/nQf  E/  =  Ei/EQ  and  I\  =  In  ,  where 

n  and  E  are  the  reservoir  values  of  the  mixture  mole-mass  ratio 
o  o 

and  molar  internal  energy  respectively,  nondimensional  forms  of 
equations  (15) ,  (16)  and  (17)  are  obtained  as 
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where,  for  the  atomic  species 


V  -V  ‘  p'tkF(1)<V2)2+  kF(!1(V"'jH 

o  *  i 

-  (L-^-nP-)  (P,)2[kR(1)  (nA*  ) 2  (nA'  )  +  kR(2)  (nA,)3J 


(21) 


Note  that  the  appearance  of  yq  in  equation  (IB)  is  strictly  due 

1/2 

to  the  arbitrary  choice  of  aQ  =  ^YoRoTo^  as  a  nondimensionalizing 

quantity;  yq  is  the  ratio  of  the  frozen  specific  heats  evaluated 

at  reservoir  conditions  and  has  no  other  significance  in  the  present 

formulation  for  a  nonequilibrium  chemically  reacting  gas. 

Equations  (18)  through  (21)  directly  yield  (|2.)  for  use  in 

2  2 

equation  (1)  where  g  =  V,  <j>  and  I\;  (  — and  are  obtained 

9 1  t 

by  differentiation  of  equations  (18)  through  (20)  with  respect 
to  t'  and  x'  respectively. 
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dwA,  3k 

These  differentiations  introduce  new  terms  such  as  t~t  ,  — r 

01  O 

3R  <CV£  3E^ 

g-jrr  /  ■  i  ,  /  as  well  as  the  corresponding  x ' -derivatives . 

Straightforward  but  somewhat  lengthy  expressions  for  these  terms 
are  derived  as  follows. 
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In  the  above,  the  electronic  energy  of  each  species  has  been  assumed 
negligibly  small  in  comparison  to  the  translational,  rotational  and 
vibrational  energy.  This  is  a  reasonable  assumption  for  the 
temperature  range  considered  in  the  present  results. 

Following  Reference  2,  the  reverse  reaction  rate  constants 
k^D  and  k^2^  are  assumed  to  vary  as 


kR1)  =  kR1)  (T,)"S 
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Note  that 


(2) 

R 

TP~ 


Sk 


=  0  when  s  =  1.  Also,  the  forward  reaction  rate 


constants  and  their  derivatives  are  obtained  from  the  relation 


k_/k_  =  K  (T) ,  where  K  (T)  is  the  equilibrium  constant. 

r  K  G  6 

Initial  conditions  at  t  =  0  are  assumed  for  p,  u,  and  T.  The 
initial  distribution  for  nA  in  the  subsonic  and  throat  regions  is 
obtained  from  a  chemical  equilibrium  calculation  using  the  local 
assumed  initial  distributions  of  p  and  T,  whereas  at  some  arbitrary 
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point  downstream  of  the  throat,  nA  is  set  equal  to  a  constant  for 
the  remainder  of  the  nozzle.  Consequently,  in  the  same  vein  as 
the  vibrational  nonequilibrium  case,  the  initial  nA  distribution 
is  at  least  qualitatively  similar  to  the  actual  steady-state 
nonequilibrium  distribution. 

Comments  on  the  Numerical  Solution 

Several  comments  are  in  order  regarding  the  numerical  behavior 
of  the  present  time-dependent  method.  First,  in  equations  (7) , 

(8),  (9),  (11),  (13),  (18),  (19)  and  (20),  the  x ' -derivatives  are 
computed  from  central  differences 


3g  _  g (x' +  Ax ' ) -g (x* -Ax ' ) 

7: sr  =  - m? - 

9  2g  _  g (x ' +Ax 1 ) -2g (x 1 ) +g (x 1 -Ax 1 ) 
3x ' 2  (Ax’)2 


evaluated  from  the  known  flow  field  at  time  t'.  Second,  the  initial 
conditions  appear  to  have  no  meaningful  effect  upon  either  the 
stability  or  convergence  of  the  technique;  results  have  been 
obtained  for  identical  nozzles  and  reservoir  conditions  but  with 
widely  different  initial  conditions,  and  in  each  case  the  same 
steady-state  solution  was  obtained.  Also,  second-order  accuracy 
as  embodied  in  equation  (1)  appears  to  be  absolutely  necessary 
for  stability.  Several  attempts  at  first-order  accuracy  (using 
the  first  two  terms  in  equation  (1) )  met  with  instabilities  for 

both  the  calorically  perfect  gas  and  the  nonequilibrium  cases. 

Hence,  the  second-order  term  involving  — — %  in  equation  (1)  is 

3 1^ 

apparently  necessary  for  convergence. 
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Regarding  the  vibrational  nonequilibrium  case,  the  nonequilibrium 
expansion  of  a  mixture  of  gases  may  result  in  complicated  expressions 
for  t  as  compared  to  the  relatively  simple  variation  assumed  above. 
Consequently,  analytical  expressions  for  and  as  presented 
above  may  be  extremely  tedious,  if  possible  at  all.  However, 
experience  with  the  present  analysis  has  shown  that  these  derivatives 
can  be  approximated  directly  by  the  finite-difference  expressions 

,3t\  _  (T')t'“(T')(t'-itl) 

(Ttr)t,  aF 


_  ,T ' (x ' +AX ' ) -T ' (x ' —Ax ' ) , 
(3xr)t,  =  1  2AX1  1 


evaluated  from  the  known  flow  field  at  times  t'  and  (t'-At'). 

The  use  of  these  finite-difference  expressions  instead  of  analytical 

1 

forms  for  and  ^p-  results  in  less  than  a  0.2  percent  change  in 
the  final  steady-state  values  for  e^^  through  the  nozzle,  and 
no  noticeable  change  at  all  in  the  other  nozzle  flow  variables. 

The  time  increment  At  appearing  in  equation  (1)  is  chosen 
at  the  end  of  each  time  step  to  satisfy  both  the  Courant-Friedrichs- 
Lewy  stability  criterion,  where 


At'  S 


Ax 1 
u '  +  a' 


(22) 


and  a  criterion  geared  to  the  speed  of  the  nonequilibrium  relaxation 
process 

At'  <  BY  (23) 
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9w ' 

where  f  =  t'  for  vibrational  nonequilibrium,  f  =  r  for  chemical 

15  a 

nonequilibrium,  and  B  is  a  proportionality  constant,  found 

empirically  to  be  approximately  0.2  for  the  present  investigation. 

Both  criteria  must  be  satisfied  simultaneously  because  two 

"relaxation"  processes  prevail  in  the  time-dependent  flow  as  it 

progresses  towards  the  steady-state;  a  fluid  dynamic  relaxation 

(eq.  (22))  involving  the  propagation  and  interaction  of  compression 

and  expansion  waves  throughout  the  nozzle,  and  an  internal  physical- 

chemical  relaxation  (eq.  (23))  involving  the  finite  rate  processes 

of  vibrational  energy  and  chemical  composition  changes. 

Because  of  the  central  differences  employed  in  the  present 
analysis,  the  flow-field  variables  at  the  first  and  last  points 
in  the  grid  (see  Figure  1)  can  not  be  directly  obtained  from 
equation  (1).  Instead,  at  the  last  point  (nozzle  exit)  all  the 
flow  properties  (in  terms  of  Z,  V,  <p ,  etc.)  are  simply  obtained 
from  linear  extrapolation  from  the  two  previous  internal  points. 

The  first  grid  point  (nozzle  inlet)  is  considered  to  be  effectively 
in  the  reservoir,  i.e.,  A/A*  at  this  point  is  a  large  number, 
usually  greater  than  ten.  Hence,  values  of  p,  T  and  p  at  this 
point  are  assumed  to  be  reservoir  values  and  held  fixed,  invariant 
with  time.  The  flow  velocity  (in  terms  of  V)  at  this  point  is 
allowed  to  vary  with  time,  and  is  found  by  linear  extrapolation 
from  the  second  and  third  internal  points.  At  large  values  of 
time,  the  velocity  at  the  inlet  approaches  its  proper  steady-state 
value  (  a  very  small  but  finite  number) ,  compatible  with  the  proper 
steady-state  mass  flow  through  the  nozzle. 
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As  a  general  comment,  the  philosophy  of  the  present  analysis 
is  similar  to  the  time-dependent  solutions  of  the  supersonic  blunt- 
body  problem  presented  in  References  16  and  17.  Because  the 
unsteady  gas-dynamic  conservation  equations  are  hyperbolic  with 
respect  to  time,  the  time-dependent  method  is  admirably  suited 
for  the  unified  analysis  of  mixed  subsonic-supersonic  flow  fields. 
Indeed,  the  expansion  through  a  convergent-divergent  nozzle  is  such 
a  case. 

RESULTS 

Using  the  present  time-dependent  analysis,  numerical  results 
have  been  obtained  on  an  IBM  7090  digital  computer  for  the  three 
cases  of  a  calorically  perfect  gas,  vibrational  nonequilibrium 
expansion  of  diatomic  nitrogen,  and  the  chemical  nonequilibrium 
expansion  of  partially  dissociated  oxygen.  Again,  emphasis  is 
made  that  the  desired  results  are  the  steady-state  flow  fields, 
and  that  the  present  time-dependent  analysis  is  an  advantageous 
means  to  that  end.  As  schematically  shown  in  Figure  1,  the 
computations  are  made  at  equally  spaced  grid  points  along  the  nozzle. 
As  few  as  30  grid  points  are  employed  in  a  calculation;  this  is 
a  convenient  number  and  does  not  represent  a  required  minimum  for 
stability.  Indeed,  a  virtue  of  the  present  time-dependent  analysis 
is  that  relatively  large  spacings  between  grid  points  can  be 
employed  throughout  the  nozzle,  including  the  near  equilibrium 
subsonic  region  near  the  reservoir.  In  addition,  using  the  present 
time-dependent  analysis  for  the  case  of  vibrational  nonequilibrium, 
several  numerical  experiments  have  been  performed  to  examine  the 
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effect  of  the  shape  (area  distribution)  of  the  subsonic  nozzle 
section  on  the  resulting  frozen  vibrational  energy  at  the  nozzle 
exit.  For  some  practical  reservoir  conditions,  the  results  show 
a  noticeable  effect,  as  will  be  discussed  in  a  subsequent  section. 
Calorically  Perfect  Gas 

In  order  to  investigate  the  purely  fluid  dynamic  behavior  of 
the  present  time-dependent  analysis,  results  have  been  obtained 
for  a  calorically  perfect  gas  with  y  =  1.4;  some  of  these  results 
are  illustrated  in  Figures  2  to  4.  Starting  with  an  initial 
linear  distribution,  the  transient  profiles  of  T'  through  the 
nozzle  are  shown  at  various  time  steps  in  Figure  2.  Two  important 
points  are  noted  from  Figure  2:  (1)  At  early  values  of  time,  the 

profiles  relax  very  rapidly  to  a  steady-state  distribution,  and 
(2)  the  resulting  steady-state  distribution  shows  excellent  agree¬ 
ment  with  known  steady-state  results  obtained  from  Reference  18. 
These  results  are  complemented  by  the  transient  profiles  of  local 
mass  flow,  puA,  through  the  nozzle,  as  shown  in  Figure  3.  The 
somewhat  wavy  initial  distribution  for  puA  is  a  consequence  of  the 
arbitrarily  assumed  initial  distributions  for  p  and  u.  Figure  3 
markedly  illustrates  that  the  transient  solution  rapidly  proceeds 
to  the  proper  steady  mass  flow,  puA  =  constant.  This  aspect  of 
the  present  analysis  is  an  important  virtue  when  applied  to 
nonequilibrium  flows,  as  described  in  the  following  sections. 
Excellent  agreement  for  the  steady-state  distributions  of  p  and  u 
is  also  obtained,  as  shown  in  Figure  4. 
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Vibrational  Nonequilibrium 

Figures  5  to  9  illustrate  results  obtained  with  the  present 

time-dependent  analysis  for  the  vibrational  nonequilibrium 

expansion  of  N2*  For  these  results,  the  characteristic  vibrational 

1/3 

relaxation  time  was  assumed  to  vary  as  ip  =  a[exp(-bT  )], 

following  Reference  11.  With  one  exception,  the  present  results 

are  obtained  with  a  =  14.7  and  b  =  0.915,  where  p  is  in  atm,  T  in 

°K,  and  t  in  seconds.  These  values  for  a  and  b  are  fitted  to  the 

vibrational  relaxation  times  used  in  References  8  and  13,  which 

19 

reflect  the  shock -tube  data  of  Blackman.  Of  course,  recent 

experiments  have  shown  that  vibrational  relaxation  times  measured 

behind  normal  shock  waves  and  subsequently  used  in  equation  (12) 

1  20 

are  not  appropriate  for  expanding  flows,  '  however,  this  does 

not  detract  from  the  present  results,  which  are  intended  to 

illustrate  the  present  time-dependent  analysis  and  which  are 

compared  with  earlier  analyses  using  Blackman's  shock-tube  data. 

To  illustrate  the  time-dependent  behavior  of  the  vibrational 

nonequilibrium  solution.  Figure  5  shows  the  transient  e^^  profiles 

at  various  time  steps;  as  in  the  previous  case,  a  rapid  approach 

to  the  steady-state  distribution  is  observed.  This  steady-state 

distribution  agrees  with  the  results  of  a  recent  steady-flow 
.  13 

analysis  by  Wilson  et  al,  as  shown  in  terms  of  T  ^  in  Figure  6. 
Both  the  present  time-dependent  analysis  and  the  steady-flow 
analysis  of  Reference  13  include  nonequilibrium  effects  upstream 
of  the  throat.  The  comparisons  in  Figure  6  are  made  for  several 
reservoir  temperatures,  and  show  very  good  agreement  between  the 
time-dependent  and  steady-flow  approaches?  however,  emphasis  is 
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again  made  that  the  present  time-dependent  results  were  obtained 
in  a  straightforward  manner,  whereas  the  steady-flow  analysis  of 
Wilson  et  al  requires  an  iterative  procedure  as  well  as  a  forward 
integration  matched  to  an  expansion  about  a  singular  point  slightly 
downstream  of  the  throat. 

For  some  practical  reservoir  conditions,  nonequilibrium  effects 

upstream  of  the  throat  are  noticeable  and  have  a  subsequent  influence 

on  the  frozen  vibrational  energy  at  the  nozzle  exit;  these  effects 

are  shown  in  Figures  7  to  9.  Figure  7  compares  results  for 

distributions  of  T  and  T,  obtained  with  the  present  time- 

vib  trans  ^ 

dependent  analysis,  which  includes  nonequilibrium  effects  upstream 
of  the  throat,  with  the  results  obtained  with  a  steady-flow  analysis 

7 

by  Harris  and  Albacete,  which  assumes  local  thermodynamic 

equilibrium  from  the  reservoir  to  the  throat.  For  this  comparison, 

1/3 

the  constants  in  xp  =  a[exp(-bT  )]  are  a  =  0.077  and  b  =  0.576, 

.  21 
which  are  fitted  to  the  quantum  mechanical  calculations  of  Widom. 

Two  important  points  are  obtained  from  the  comparison  shown  in 
Figure  7;  (1)  the  present  results  show  a  noticeable  departure 

from  equilibrium  flow  upstream  of  the  throat  as  reflected  by  the 
difference  between  Ttrans  and  Tv^  and  (2)  these  upstream  non¬ 
equilibrium  effects  result  in  a  higher  frozen  T  ^  at  the  nozzle 
exit  in  comparison  to  the  results  of  Reference  7.  For  the 
reservoir  conditions  shown  in  Figure  7,  the  vibrational  energy  is 
a  relatively  small  fraction  of  the  internal  energy  of  the  gas. 
lienee,  the  differences  in  between  the  two  analyses  are  not 

paralleled  by  similar  differences  in  Ttrans»  indeed,  good  agreement 
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v 


is  obtained  for  T  .  In  addition,  good  agreement  is  also 

tr  sns 

obtained  for  p  and  u  as  shown  in  Figure  8.  The  upstream 
nonequilibrium  effects  are  accentuated  at  lower  pressures,  as 
markedly  shown  in  Figure  9.  here ,  the  present  results  are  compared 

g 

with  the  steady-state  analysis  of  Erickson,  which  like  that  of 
Reference  7  assumes  equilibrium  flow  to  the  throat.  The  results, 
which  are  given  in  terms  of  local  nonequilibrium  vibrational  energy 
non-dimensionalized  by  the  reservoir  equilibrium  value,  again 
illustrate  the  effect  of  nonequilibrium  conditions  upstream  of  the 
throat,  particularly  for  the  lower  value  of  pQ.  Consequently, 

Figures  7  and  9  demonstrate  the  importance  of  treating  nonequi¬ 
librium  effects  upstream  of  the  throat  when  solving  the  nonequi¬ 
librium  nozzle  flow  problem.  In  contrast  to  a  steady-flow 
approach,  the  present  time-dependent  analysis  handles  such  effects 
in  a  straightforward  manner. 

Influence  of  Subsonic  Area  Distribution 

The  fact  that  significant  nonequilibrium  conditions  can 
prevail  upstream  of  the  throat  prompts  the  following  question: 

What  is  the  effect  of  the  shape  (area  distribution)  of  the  subsonic 
section  on  the  resulting  frozen  vibrational  energy  at  the  nozzle 
exit?  In  order  to  examine  this  question,  several  numerical 
experiments  hcve  been  performed  using  the  present  time-dependent 
technique.  Figures  10  and  11  reflect  some  results  from  these 
experiments  in  which  vibrational  temperature  distributions  are 
compared  for  two  nozzles  with  different  subsonic  sections  but 
identical  supersonic  sections.  These  nozzles  are  shown  schematically 
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in  Figure  10;  the  composite  nozzle  exhibits  a  greater  subsonic  area 
change  near  the  throat  than  the  purely  hyperbolic  nozzle.  Figure  10 
shows  only  a  small  difference  in  frozen  T  ^  at  the  nozzle  exit 
for  a  reservoir  pressure  of  80  atm.  However,  Figure  11  shows  that 
at  a  lower  pressure  of  10  atm  a  considerable  difference  in  frozen 
vibrational  temperature  is  obtained  and  this  difference  appears 
to  be  relatively  insensitive  to  the  reservoir  temperature. 
Consequently,  for  practical  reservoir  conditions,  some  tailoring 
of  the  frozen  vibrational  energy  at  the  nozzle  exit  can  be 
accomplished  by  proper  design  of  the  subsonic  section. 


where  =  02  and  M2  =  0.  The  reaction  rate  constants  and  nozzle 

.  1  2 
shape  are  taken  identical  to  those  of  Hall  and  Russo,  '  where 

kRX)  =  kR1)  < VT)  '  kR2>/kR1)  =  35  T/TD'  TD  =  D/k'  D  is  the 
o 

dissociation  energy  per  molecule,  and  kRVkR3  =  Ke<  For  the 
present  results,  the  equilibrium  constant,  Kg,  is  obtained  from 
Wray33  as  K  =  (1.2  x  103)T  ^3  exp(-118,000/RT)  ,  where  K  is  in 

C  6 

moles/cm^  and  T  in  °K.  Also  in  the  present  results  the 

vibrational  energy  is  assumed  to  be  in  the  local  thermodynamic 

equilibrium  at  the  local  gas  temperature. 
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Figure  12  illustrates  the  transient  profiles  of  atom  mass 
fraction  a  at  various  time  steps,  where  a  =  p^/p  =  As  in 

the  previous  cases,  a  continuous  approach  to  the  steady  state  is 
observed.  However,  in  comparison  with  Figures  2  and  5  from  the 
previous  cases,  Figure  12  shows  that  the  steady  state  for  the 
chemical  nonequilibrium  case  is  obtained  after  considerably  more 
time  steps.  This  behavior  is  due  to  the  small  values  of  At 
dictated  by  equation  (23)  for  the  present  chemical  nonequilibrium 
conditions;  this  contrasts  to  the  previous  results  for  the 


calorically  perfect  gas  and  vibrational  nonequilibrium  cases, 
where  the  minimum  At  came  from  equation  (22)  and  where  the  values 
of  At  were  an  order  of  magnitude  larger  than  the  present  chemical 
nonequilibrium  case.  Also  in  Figure  12,  the  resulting  steady- 
state  distribution  is  compared  to  the  results  obtained  from  the 
steady-flow  analysis  by  Hall  and  Russo;  very  good  agreement  is 
obtained.  (Note  that  a  slight  discrepancy  occurs  with  the 
equilibrium  reservoir  mass  fraction,  aQ,  which  can  be  attributed 
to  slight  differences  in  Kg  between  the  two  sets  of  results.) 

Whereas  the  results  shown  in  Figure  12  are  computed  for 
pQ  =  9.4  atm.  Figures  13  to  17  show  results  obtained  for 
pQ  =  82  atm,  where  equilibrium  conditions  are  more  closely 
approached,  hence  providing  a  more  stringent  test  for  the  stability 
of  a  nonequilibrium  flow  analysis.  Figure  13  shows  the  variation 
of  a  at  the  nozzle  exit  as  a  function  of  time;  the  steady-state 


aexit  aPProac^e^  at  large  values  of  time.  Figures  14  to  17 


show  results  obtained  with  the  present  time-dependent  analysis 
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for  the  steady-state  a,  T,  p  and  p  distributions,  respectively 
(plotted  versus  nozzle  area  ratio) .  These  results  are  compared 
with  the  results  of  Hall  and  Russo,  obtained  from  a  steady-flow 
analysis.  Very  good  agreement  is  obtained;  however,  emphasis 
is  again  made  that  the  present  results  are  obtained  in  a 
straightforward  manner,  whereas  the  steady-state  analysis  of 
Reference  2  requires  an  iterative  procedure  due  to  the  unknown 
critical  mass  flow,  and  employs  an  asymptotic  series  solution 
for  the  flow  at  large  areas  upstream  of  the  throat  for  starting 
the  forward  integration. 

CONCLUSIONS 

In  summary,  the  following  conclusions  are  made:  (1)  A  new, 
alternative  solution  for  quasi-one-dimensional  nonequilibrium 
nozzle  flows  is  presented.  A  virtue  of  this  new  technique  is 
its  simplicity,  which  prevails  from  its  initial  physical 
formulation  to  the  successful  receipt  of  numerical  results. 
Consequently,  the  time-dependent  technique  appears  to  warrant 
consideration  for  future  applications  in  nonequilibrium  nozzle 
flows.  (2)  Because  in  some  practical  cases  nonequilibrium 
conditions  prevail  upstream  of  the  throat,  the  shape  of  the 
subsonic  section  can  have  a  noticeable  effect  on  the  frozen 
vibrational  energy  at  the  nozzle  exit. 

Parenthetically,  it  should  be  noted  that  the  present  time- 
dependent  results  always  converged  to  the  supersonic  solution 
downstream  of  the  throat;  a  purely  subsonic  solution  was  never 

observed.  This  situation  was  due  to  the  choice  of  initial  values 
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which  were  in  qualitative  agreement  with  a  supersonic  solution 
(i.e.,  low  exit  values  of  density  and  temperature,  and  high 
values  of  velocity) . 
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FIG.  3  TRANSIENT  AND  FINAL  STEADY -ST ATE  MASS-FLOW  DISTRIBUTIONS 
FOR  A  CALORICALLY  PERFECT  GAS  OBTAINED  FROM  THE  PRESENT 
TIME-DEPENDENT  ANALYSIS;  Y=  1.4 
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FIG.  4  STEADY -STATE  PRESSURE  AND  VELOCITY  DISTRIBUTIONS 
FOR  A  CALORICALLY  PERFECT  GAS;  COMPARISONS  OF 
THE  PRESENT  TIME  -  DEPENDENT  ANALYSIS  WITH  RESULTS 
OBTAINED  FROM  REFERENCE  18. 


TRANSIENT  AND  FINAL  STEADY-STATE  eV|B  DISTRIBUTIONS 
FOR  THE  NONEQUILIBRIUM  EXPANSION  OF  N2  OBTAINED 
FROM  THE  PRESENT  TIME -DEPENDENT  ANALYSIS 


VIBRATIONAL  TEMPERATURE,  TV|B,  (°K  xIO' 


NOLTR  69-52 


FIG.  6  STEADY -STATE  Tv|6  DISTRIBUTIONS  FOR  THE  NONEQUILIBRIUM 
EXPANSION  OF  N2;  COMPARISON  OF  THE  PRESENT  TIME- 
DEPENDENT  ANALYSIS  WITH  THE  STEADY-FLOW  ANALYSIS  OF 
REFERENCE  13. 
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FIG.  7  STEADY-STATE  TV|B  AND  TjranS  DISTRIBUTIONS  FOR  THE  NONEQUILIBRIUM 
EXPANSION  OF  N2;  COMPARISON  OF  THE  PRESENT  TIME-DEPENDENT  ANALY¬ 
SIS  WITH  THE  STEADY -FLOW  ANALYSIS  OF  REFERENCE  7. 


STEADY  FLOW  ANALYSIS  OF 
ERICKSON  (REF.  8) 
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FIG.  9  STEADY -ST  ATE  eV(B  DISTRIBUTIONS  FOR  THE  NONEQUILIBRIUM  EXPANSION 
OF  N2;  COMPARISON  OF  THE  PRESENT  TIME-DEPENDENT  ANALYSIS  WITH 
THE  STEADY -FLOW  ANALYSIS  OF  REFERENCES. 


FIG.  10  COMPARISON  OF  TV(B  DISTRIBUTIONS  FOR  TWO  NOZZLES  WITH 
DIFFERENT  SUBSONIC  SECTIONS  BUT  IDENTICAL  SUPERSONIC 
SECTIONS;  P0  80  atm,  NONEQUILIBRIUM  EXPANSION  OF  N2 


VIBRATIONAL  TEMPERATURE,  (°K)  x  10" 
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FIG.  11  COMPARISON  OF  TvjB  DISTRIBUTIONS  FOR  TWO  NOZZLES  WITH 
DIFFERENT  SUBSONIC  SECTIONS  BUT  IDENTIC/  L  SUPERSONIC 
SECTIONS;  p0  10  atm,  NONEQUILIBRIUM  EXPANSION  OF  N2. 
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FIG.  12  TRANSIENT  AND  FINAL  STEADY-STATE  ATOM  MASS  FRACTION  DISTRIBUTIONS 
FOR  THE  NONEQUILIBRIUM  EXPANSION  OF  DISSOCIATING  OXYGEN  OB¬ 
TAINED  FROM  THE  PRESENT  TIME-DEPENDENT  METHOD;  THE  STEADY-STATE 
DISTRIBUTION  IS  COMPARED  WITH  THE  STEADY-FLOW  ANALYSIS  OF  REFERENCE  2. 
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FIG.  14  STEADY-STATE  ATOM  MASS  FRACTION  DISTRIBUTION  FOR 
THE  NONEQUILIBRIUM  EXPANSION  OF  DISSOCIATING 
OXYGEN;  COMPARISON  OF  THE  PRESENT  TIME-DEPENDENT 
ANALYSIS  WIIH  THE  STEADY  -FLOW  ANALYSIS  OF  REF.  2 
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FIG,  15  STEADY -STATE  TEMPERATURE  DISTRIBUTION  FOR  THE  NON¬ 
EQUILIBRIUM  EXPANSION  OF  DISSOCIATING  OXYGEN; 
COMPARISON  OF  THE  PRESENT  TIME-DEPENDENT  ANALYSIS 
WITH  THE  STEADY-FLOW  ANALYSIS  OF  REF.  2 


NON-DIMENSIONAL  DENSITY:  ppQ 


NON-DIMENSIONAL  PRESSURE:  p/p 
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FIG.  17  STEADY-STATE  PRESSURE  DISTRIBUTION  FOR  THE  NON¬ 
EQUILIBRIUM  EXPANSION  OF  DISSOCIATING  OXYGEN; 
COMPARISON  OF  THE  PRESENT  TIME-DEPENDENT  ANALY¬ 
SIS  WITH  THE  STEADY  FLOW  ANALYSIS  OF  REF.  2 
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A  a  b  TRAC  r 


SPONSORING  MILITARY  A  C  T|  /  I  Tv 


A  new  technique  is  presented  for  the  numerical  solution  of 
quasi-one-dimensional ,  vibrational  and  chemical  nonequilibrium 
nozzle  flows  including  nonequilibrium  conditions  both  upstream 
and  downstream  of  the  throat.  This  new  technique  is  a  time- 
dependent  analysis  which  entails  the  explicite  finite-difference 
solution  of  the  quasi-one-dimensional  unsteady  flow  equations  in 
steps  of  time,  starting  with  assumed  initial  distributions 
throughout  the  nozzle.  The  steady-state  solution  is  approached 
at  large  values  of  time.  A  virtue  of  the  present  time-dependent 
analysis  is  its  simplicity,  which  prevails  from  its  initial  physical 
formulation  to  the  successful  receipt  of  numerical  results.  Also, 
the  present  solution  yields  the  transient  as  well  as  the  steady- 
state  nonequilibrium  nozzle  flows.  To  exemplify  the  present  analysis, 
results  are  given  for  several  cases  of  vibrational  and  chemical 
nonequilibrium  expansions  through  nozzles. 


DD  F0RM  1473 


(PAM 


UNCLASSIFIED 

S<-<  'ini'.-  iilion 


Nonequilibrium  flow 
Time-dependent  analysis 
Chemical  nozzle  flows 
Vibrational  nozzle  flow 
Nozzle  expansions 


